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Cramer-Rao Bound for Localization with A 
Priori Knowledge on Biased Range 
Measurements 

Tao Wang 
Abstract 

This paper derives a general expression for the Cramer-Rao bound (CRB) of wireless localization 
algorithms using range measurements subject to bias corruption. Specifically, the a priori knowledge 
about which range measurements are biased, and the probability density functions (PDF) of the biases 
are assumed to be available. For each range measurement, the error due to estimating the time-of-arrival 
of the detected signal is modeled as a Gaussian distributed random variable with zero mean and known 
variance. In general, the derived CRB expression can be evaluated numerically. An approximate CRB 
expression is also derived when the bias PDF is very informative. Using these CRB expressions, we 
study the impact of the bias distribution on the mean square error (MSE) bound corresponding to the 
CRB. The analysis is corroborated by numerical experiments. 

Index Terms 

Cramer-Rao bound, localization, biased range measurements, wireless networks. 

I. Introduction 

Wireless localization systems have been attracting intensive research interest in both academia and 
industry lately. For indoor or dense urban environments, Global Positioning Systems do not function well 
for geolocation purposes. Instead, a few beacons at fixed and known positions are exploited for geolocation 
using measurements related to signal power decay, bearing, difference of range, and range between a 
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target and the beacons JT]. Thanks to the superior multipath resolution and penetration capability of ultra- 
wideband (UWB), a range measurement can be obtained by UWB pulsed signals with high accuracy GJ. 
Therefore, range based localization has become a promising option for short-range wireless networks. 

Much work on ranging algorithms based on time of arrival (TOA) measurements has been reported 
El . |[3l . Specifically, they are usually based on detecting the first arriving signal from a transmitter 
and estimating its time of arrival (TOA) at a receiver. Note that the clocks of the transmitter and the 
receiver need to synchronized, and the transmission time at the transmitter should be sent to the receiver 
for producing the range measurements. In this paper, we consider the geolocation algorithms using 
range measurements generated between the target and multiple beacons in the aforementioned way. It 
is important to note that the considered algorithms are suitable for cooperative (blue-force) geolocation 
applications, in that they require the cooperation of the target and the beacons in order to synchronize 
their clocks for generating range measurements. 

In general, a range measurement error can be modeled as the sum of two terms. The first one is due to 
the error of estimating the TOA of the detected signal, while the second one corresponds to the difference 
between the path length traveled by the detected signal and the transmitter-receiver distance. For instance, 
the second term is nonzero when the detected signal does not come from the line of sight propagation 
path. The first term is usually modeled as a zero mean random variable with a Gaussian distribution in 
the literature EJ-JSJ. As a consequence, the second term is equal to the bias of the range measurement. 

It is well known that the Cramer-Rao bound (CRB) sets a fundamental lower limit to the co variance 
of any unbiased estimator for a vector parameter (please refer to pages 63 to 72 in Q). Specifically, 
the mean square error (MSE) of any unbiased location estimator is lower bounded by the corresponding 
CRB's trace, referred to as the MSE bound hereafter. The CRB and the MSE bound are used extensively 
to evaluate particular localization algorithms, as well as guide the geolocation system design fulfilling 
certain accuracy requirements |@], ||8], |0. When some of the range measurements are biased, the CRB 
has been derived in ifTOl by regarding each bias as a deterministic nuisance parameter, and it was shown 
in O that using those bias-corrupted measurements does not improve the CRB. 

Recently, much work has been reported on modeling the bias as a random variable with a particular 
probability density function (PDF). Specifically, this PDF is obtained empirically by extensive experiments 
ifTTIl . lfT2l . or theoretically by scattering models |fl"3l . In case the PDF of the bias is known a priori, a 
Cramer-Rao-like bound has been derived in [ 14j for the joint estimation of the deterministic target location 
and the random bias. This bound is referred to as the hybrid CRB in |fl"5ll , lfT6ll , since the parameters 
to be estimated consist of both deterministic and random parameters. It was shown in |5j that using the 
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bias-corrupted measurements improves the MSE bound associated with the hybrid CRB. 

Specifically, the hybrid CRB was derived from the joint PDF of the distance measurements and the 
random bias. Compared with the hybrid CRB, the CRB derived from the marginal PDF of range 
measurements is much tighter and can be asymptotically attained by a maximum likelihood (ML) 
estimator 0, COD. When the bias PDF is approximated by a piecewise constant function, an expression 
that can be numerically evaluated was derived for the CRB in |6]. Numerical results showed that for 
uniformly distributed bias, the presence of the bias degrades the MSE bound and using the bias-corrupted 
measurements improves the MSE bound. 

Compared with the above existing work, this paper contains the following contributions: 

• We derive a general expression that can be numerically evaluated for the CRB. When the prior bias 
PDF is very informative, an approximate CRB expression is derived. 

• Based on these expressions, the impact of the bias distribution on the MSE bound corresponding to 
the trace of the CRB is studied analytically. 

The rest of this paper is organized as follows. Section II describes a typical 2D localization system 
and related models for range measurements. In Section III, a general CRB expression is derived for this 
system. After that, an approximate CRB expression is derived in Section IV when the prior bias PDF is 
very informative. Based on these expressions, the impact of the bias distribution on the MSE bound is 
studied analytically in Section V. We will show some numerical results in Section VI, and complete this 
paper by some conclusions in Section VII. 

Notations'. Upper (lower) boldface letters denote matrices (column vectors), and [■] represents the 
transpose operator. V x (/(x)) stands for a column vector which is the gradient of the function /(x) with 
respect to x, and V^.(/(x)) represents a matrix which is the Hessian of the function /(x) with respect 
to x. 

II. SYSTEM SETUP AND RANGING MODELS 

We consider a typical 2D geolocation system equipped with M beacons (M > 3). In the following 
sections, although only the CRB for this 2D system is derived, the CRB for a 3D system can be derived in 
the same way and will be given as well. The coordinate of beacon m is p m = [x m , y m } T (m = !,■■■ ,M), 
and a target is located at u = [x u ,y u ] T . 

Suppose independent range measurements have been produced between a target and beacons, and the 
one between beacon m and the target is denoted by r m . Without loss of generality, we assume the first L 
range measurements are biased, while the other measurements are all unbiased. The measurement error 
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(h, 



m = 1, ■ ■ ■ , L 
m = L + !,-■■ ,M, 



(1) 



related to r m is modeled as: 

where d m represents the distance between the target and beacon m, and e m is due to the TOA estimation 
error of the detected signal, and b m corresponds to the difference of the path length traveled by the 
detected signal and the distance from beacon m to the target. We assume b m has the a priori known 
PDF p(b m ) with its mean and variance denoted by b m and K m , respectively. In addition, we assume b m 
is independent of e m , and e rn is Gaussian distributed with zero mean. Besides, we denote its variance 
as a m and assume it is known a priori. The motivation behind this assumption is twofold. One is that 
it simplifies the mathematical derivation. The other is that the impact of the bias distribution on the 
localization CRB can be studied when the TOA estimation has a guaranteed accuracy at a prescribed 
level. 

We stack all range measurements into the column vector r = [r±, ■ ■ ■ , tm] T ■ The log-likelihood function 
of u is denoted as A = ln(p(r;u)) where p(r; u) is the PDF of r. Since all range measurements are 
independent of each other, A = J2m=i ^-m, where A m = ln(p(r m ; u)) is the log-likelihood function of 
r m given u, and p(r m ; u) is the PDF of r m . 

For an unbiased range measurement r m (m = L + 1, • • ■ , M), p(r m ; u) can be formulated as: 



p(r m ;u) 



1 



: exp 



2al 



(2) 



For a biased range measurement r m (m = 1, • • • ,L), we define T m = ln(p(r m \b m ; u)) as the joint 
log-likelihood function of r m given u and b m , where p(r m \b m ; u) is the PDF of r m conditional on b m 
and u. p(r m \b m ; u) and p(r m ; u) can be expressed respectively as: 



p{r m \b m ;u) 

p{r m ;u) 



1 



exp 



\/ 27rc7 n 

p{r m \b m ]\i)p(b m )&b 



(3) 
(4) 



III. Derivations of a general crb expression 

The CRB is the inverse of the Fisher information matrix (FIM) F, which is expressed by: 

F = E r;u [(V u A)(V u A) T ] (5) 

where E r;u [■] denotes the expectation operator with respect to p(r; u). Since all distance measurements 
are independent of each other, F = ^ . F m where F m is the FIM related to r m Q: 



Er m ;u [(V u A m )(V u A„ 



(6) 
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where E r . m;u [•] denotes the expectation operator with respect to p(r m ; u). 
When m = L + 1, • • ■ , M, we can show that: 



-<t> m 4>m 



(V) 



where (j> r 



When m = 1 , ■ 



represents the orientation of the target relative to beacon m. 
, L, we can first prove that: 
1 



p(r m ;u 

30 f>(Tm|&m;u)p(6 m ) 



E fe m |r m ;u [V u r m ] 



[V u p(r m |6 m ; u)] p(b m )db r , 
[V u r m ]d6 m 



°"m a m4 > r l 



(8) 



where a m = E fem | rm . u [r m — d m — b m ] and E bm | rm . u [-] denotes the expectation operator with respect 
to the posterior PDF p(b m \r m ;u). This means that a m actually stands for the posterior mean of e m 
conditional on r m . Inserting (HJ into ©, we can compute F m by: 



F m — A m m m 

where A m = a~ 4 ■ E rm;u [a^] . 

Summing Fi, • • ■ , Fjv/ and taking the inverse, we can express the CRB as: 

" L M 

ET ^ — ^ 2 T. 

\ m (j> m 4> m + 2^ °m ^m^n 



(9) 



(10) 



_m=l m=L+l 

The CRB for a 3-D localization system can be derived in the same way as given above. Most 
interestingly, we can show after simple mathematical arrangement that the CRB for the 3-D localization 
system can still be evaluated by (flOl) . when u and p m are substituted with the 3-D coordinates of the 
target and beacon m, respectively. 

It is important to note that (fTOl is a general expression for the CRB, when each biased measurement 
r m is Gaussian distributed conditional on b m . However, we have proven that V u A m = E& i rm . u [V u r m ] 
holds in general for any particular form of p{r m \b m ; u). This means that a closed-form CRB expression 
can be derived in a similar way when p{r m \b m ; u) takes other specific form. 

To illustrate the generality of (fTOl . we show that the CRB can be derived by (fTOb when b m is determinis- 
tic. Note that this CRB will be used when illustrating the effectiveness of the approximate CRB derived in 
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Section IV. In this case, p{b m \r m ; u) = S(b m -b m ), a m = r m -d m -b m , p{r m ; u) = p{r m \b m = b m ; u), 
and therefore A m = cr m 2 . This means that the contribution of r m to the CRB is the same as if this range 
measurement were not bias corrupted. This is true since when b m is deterministic and known a priori, 
its value can be subtracted from r m , which means that r m is in effect not subject to the bias corruption. 

In general, an analytical expression for A m does not exist, in such a case A m can be evaluated 
numerically, with C following from equation (10). This provides a way to study the effects of the 
geometric configuration, the TOA estimation quality, as well as the bias distribution on the CRB. To this 
end, it can easily be shown that 

/+oo 
o^(r m ;u)dr m (11) 
-oo 

and 

f +Q ° I A h M r m\b m ;u)p(b m ) 

a m = {r m -d m -b m ) — db m , (12) 

J-oo P\ r m',' 11 ) 

respectively. For every value of r m , p(r m ;u) and a rn can be first evaluated by computing numerically 
the integrals in dU and (ITZb with respect to b m Q, respectively. Then, A m can be evaluated by computing 
numerically the integral with respect to r m in (flTT ). 

Note that the CRB expression in (fTOl . as those presented in iffl-O, is derived under the assumption 
that we have the a priori knowledge about which range measurements are biased. In almost all practical 
scenarios, this a priori knowledge may not be available lTT8l - ll20l . In this case, the CRB in (flOl) is an 
optimistic lower bound, since it assumes the implicit knowledge which is not really available for the 
range measurements. We will illustrate the optimism of (TTOl when the a priori knowledge about which 
range measurements are biased is not available by numerical experiments in Section VI. Relaxing the 
assumption of knowing which measurements are biased, and generating an information reduction factor 
that scales the CRB when the probability of each measurement getting biased is known a priori, in a 
similar way as shown in |[T8l - ll20l . would be a valuable piece of future work. 

We should also note that the derived CRB is for cooperative geolocation applications, since the 
synchronization between the target and the beacons is required as mentioned in the introduction part. For 
geolocating evasive targets, range measurements would only be available with active sensing, in which 
case azimuth measurements, etc., would be available as well. It is another interesting piece of future work 
to derive the CRB for localization using azimuth measurements, etc., in addition to range measurements. 

'The detailed procedure of numerical integration can be found in pages 198 — 200 in 1171 . 
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IV. Derivation of an approximate crb expression when n m < a m 

We now derive an approximate CRB expression when b m has a small variance compared to a m (i.e. 
K m <C (T m ). First of all, we express F. m (m = 1, • • • , L) alternatively as Q: 

F m = -E rm . u [Vl A m ] (13) 

Next, we approximate p(r m |6 m ;u) by the first-order Taylor-series expansion as follows: 

p(r m \b m ; u) « / m + f' m ■ (b m - b m ) (14) 

where f m and f m are expressed respectively as: 



fm = p(r m \b m = b m ;u) 
, dp(r m \b m ;u) 

Jm 



db r . 



b m =b„ 



(15) 
(16) 



Based on (01) and (fT4l ). we can show that p(r m ;u) « / m . In addition, we make the following 
approximations: 



P(r m \bm] u) ^ /m + /m ■ (fern ~ frn 

p(r TO ;u) 



= 1 + 



= 1 + 



db r . 



b m =b, 

dm b* 



' m ""771 u m 



{b m ~ b m ) 
{bm b m ) 



(17) 



Inserting (fTTT ) into ([T2l . we can approximate a m as follows: 



l _ 1m 
cri. 



(Xm dm bm) 



Based on ([8]) and ( fT8l ), V„A m can be approximated by: 

°m 2 [(V u a m )0m + a m (Vu4)] 



I" ^ ) <i>m<t> m + am{Vid. 



Using (fTST ). (fl9l . and (U~3l . F m can be approximated by: 



T 



where 



A m = °"m 2 (l - (Km/Vm) 2 ) ■ 



(18) 



(19) 



(20) 



(21) 
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From © and <EQ}, we can see that A 



•in 



In addition, CRB can be approximated by: 



r l m 



-, -l 



(22) 



_m=l m=L+l 



Note that a better approximation of the CRB can be derived with a higher order Taylor-series expansion 
in the same way as given above, when the statistical moments of b m higher than the 2nd order are known 
a priori. Here, we only use the first order Taylor-series expansion to derive an approximation when the 
2nd order moment of b m (i.e. K m ) is known. The derived approximate CRB expression features a simple 
mathematical structure but reveals the influence of the bias distribution on the CRB clearly. 

We will illustrate the effectiveness of d22l by numerical experiments using measured bias distributions in 
Section V. In fact, we can analytically illustrate the effectiveness of (1221 ) when b rn is Gaussian distributed. 
Note that although in practice the Gaussian distribution is unlikely to be a good approximation of the 
true bias distribution, it lends the CRB to be easily simplified analytically so that the effectiveness of 
(1221 ) can be examined. In this case, r m can be equivalently modeled as r rn = d m + b rn + e' m where e' m is 
Gaussian distributed with zero mean and variance cr^ + F? m . Therefore, A m = (cr^ + k^) -1 according 
to the earlier analysis for the case when the bias is deterministic in Section III. Since n m <C a m , we can 
see that C can indeed be approximately computed with d22l . due to the fact that 



It is important to note that the MSE bound, namely the trace of C and denoted by Tr(C), reduces 
as A m increases. This is because Ylm=i increases in a positive definite sense, therefore both C and 
Tr(C) decrease, which means that the MSE bound improves, as X m increases. 

When r m (m = 1, • • • , L) are discarded, we can compute C by setting A m to 0. Since A m is no smaller 
than 0, using a bias-corrupted measurement results in the same or better MSE bound than discarding it. 



A 



•in 




(23) 



V. MSE BOUND ANALYSIS 
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An important property of A m is that A m < a m 2 . This can be justified by: 



A 



Til 



< a" 4 • E rm;u [E 6m | rm . u [(r m - d m - 

= °"m 4 ' E b m [ E r m |6 m ;u [{ r m ~ d m - b. 



b m ]f 



-2 



(24) 



m 



where Et m [•] and E Tm \ bm . u [•] stand for the expectation operator with respect to p(b m ) and p(r m |6 m ; u), 
respectively. Specifically, the second inequality is according to the Jensen's inequality, and the equality 
holds if b m is deterministic (please refer to pages 77 — 78 in [21]). This implies the presence of the 
random bias b m degrades the MSE bound. 
It is interesting to consider two special cases: 

1) The first case is when b rn has a very large variance compared to a m (i.e., K m 3> cx m ) so that for 
any fixed r m , p(b m ) is approximately constant within the nonzero support of p{r m \b m ; u) . This 
case corresponds to the scenario where the prior bias PDF is not informative, and we can find that 
A m sa 0. If b m is Gaussian distributed, this is the case when n m is sufficiently large, and A m is 
indeed close to zero. This means that discarding a bias-corrupted measurement degrades the MSE 
bound slightly when the prior bias PDF is not informative. 

2) The second case is when b m has a very small variance compared to a m (i.e., n m <C cr m ) so that for 
any fixed r m , p(r m \b m ;u) varies slowly with respect to b m within the nonzero support of p(b m ). 
This case corresponds to the scenario where the prior bias PDF is very informative. Based on the 
analysis in Section IV, we can find that A m « X' m = a m 2 (l — (K m /cr m ) 2 ) « <r~ 2 . This means that 
the presence of the random bias b m results in a slight degradation to the MSE bound, compared to 
the case in which there are no measurement biases. 



For illustration purposes, we consider a wireless localization system with four beacons located respec- 
tively at pi = [0, 10] r , p 2 = [0, -10] T , p 3 = [-10, 0] T , p 4 = [10, Of, and a target at u = [-1, -5] T , as 
shown in FigureQ] We assume n is biased while r-i, ^3, and are unbiased, and a m = 1 (m = 1, ■ ■ • ,4). 
p(b\) is assumed to have the same shape as the measured bias distribution reported in [11]. Specifically, 



VI. Numerical experiments 
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x-axis (m) 



Fig. 1. The localization system with 4 beacons considered in the numerical experiments. Note that only ri, namely the range 
measurement between u and pi, is biased. 

TABLE I 

The value of Pi, i = 0, • • • ,8. 



Po 


Pi 


P 2 


P 3 


Pa 




Pe 


Pi 


Ps 


0.12 


0.03 


0.31 


0.12 


0.24 


0.12 


0.03 





0.03 



is expressed as 

( if 6i € (-oo,fi ]U(O 9 ,oo) 

P(&i) = S p . (25) 
[ % if &i € (n 4 , fii+i] with i = 0, • • • ,8, 

where = 0.1 + £A, and Pi is given in Table U for i = 0, • • • , 8. When A = 0.1, p{b\) is shown in 

Figure [2] It can be readily derived that b\ = 0.1 + 3.49A and k\ = 1.83A, which means that k\ is a 

linear function of A. When A increases from 0.1 to 1, we have computed k±, and the results are shown 

in Figure [3] Note that all the above values related to a coordinate or length have the units of meters. 

To examine the effectiveness of the approximate CRB expression (l22l) . we have evaluated the CRB 
from ( fTOb numerically, and the approximate CRB from d22t when A increases with m / a m < 1 satisfied. 
The MSE bounds corresponding to those CRBs are shown with respect to ni/a m in Figure HJa. We can 
see that when K\ja m < 0.5, the MSE bound corresponding to the approximate CRB is very close to the 
one computed from ( fTOb . This illustrates the effectiveness of (1221) when K\ <C a m . 

To examine the analysis in Section V, we have also evaluated the CRB from (fTOb numerically when A 
increases with K\/a m > 1 satisfied, and the corresponding MSE bound is shown with respect to Ki/a m 
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0.2 0.4 0.6 0.8 1 1.2 1.4 

bias (m) 



Fig. 2. The PDF of 61, which is the bias of the range measurement between u and pi, when A = 0.1 m. 




A(m) 

Fig. 3. Ki, which is the standard deviation of the random bias 61, when A increases from 0.1 to 1 m. 

in Figure @]b. We have also computed the CRBs when r\ is discarded and when r\ is not bias corrupted, 
respectively, and shown the corresponding MSEs in both Figure |4]a and Figure @]b. It is clearly shown in 
Figure H]a that as K\ reduces, the MSE bound approaches the one corresponding to the case when r± is 
not bias corrupted. On the other hand, as k\ increases, the MSE bound approaches the one corresponding 
to the case when r\ is discarded, as shown in Figure |4]b. These observations corroborate the analysis in 
Section V. 

Note that the derived CRB is an optimistic bound in scenarios when we do not have the a priori 
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1.4r 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

1 m 



(a) when tii/a m < 1 




1.05 :,-;-„■ 

1 5 10 15 20 25 30 35 

(b) when K\ja m > 1 

Fig. 4. The MSE bounds corresponding to the CRBs computed when Ki/a m < 1 and Ki/a m > 1, respectively. 

knowledge about which range measurements are biased. To illustrate this remark explicitly, we compare 
the MSE bound corresponding to the CRB computed from (fTOT l with the MSE of two ML location 
estimators as K\ja m increases, since the ML estimator is widely used and able to achieve the corre- 
sponding CRB asymptotically. For the first ML estimator, the above mentioned a priori knowledge is not 
available, and a binary parameter s m is introduced to indicate that r m (m = 1, ■ ■ ■ , 4) is bias corrupted 
if s m = 0. Specifically, the first ML estimator produces a joint estimate of u and {s m ,m = 1, • • • ,4} 
as the maximizer of the log-likelihood function 

4 

L(u, s±, s 2 , s 3 , s 4 ) = L m (u, s m ), (26) 

m=l 
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Fig. 5. The MSEs of the ML location estimators and the MSE bound corresponding to the CRB computed from (10). Note 
that the MSE bound is the same as the corresponding one depicted in Fig. [4] 



where 



L m (u, s m ) = ln( 1 ) + In ( s m exp{- ^ m } 
V 27rc7 m V Zcr m 

+(1 - s m ) r°°p{b m ) exp{- (rm - f m ~ hm)2 }db m ) , (27) 



2fi 2 

J — oo ^ u m / 

and p(b m ) = p{b\), m = 2,3,4. For the second ML estimator, the a priori knowledge is available, i.e., 
this estimator knows a priori that only n is biased, and produces an estimate of u as the maximizer of 
L(u, si, S2, S3, S4) when si = and s m = 1, m = 2,3,4. For every A, the MSE of each estimator is 
computed by averaging the square errors of the location estimates for 1000 random realizations of range 
measurements, and the results are shown in Figure |5]with respect to K\/a m . We can see that the MSE 
of the second ML estimator is slightly above the MSE bound computed from (TTOt . while the MSE of the 
first ML estimator is much greater than that MSE bound. These observations indicate that the derived 
CRB is indeed an optimistic lower bound for localization algorithms without the a priori knowledge 
about which range measurements are biased. 

VII. Conclusion 

We have derived a general expression for the CRB of wireless localization algorithms using range 
measurements subject to bias corruption. Specifically, the knowledge about which range measurements 
are biased, and the probability density functions (PDF) of the biases are assumed to be known a priori. 
For each range measurement, the error due to estimating the TOA of the detected signal is modeled 
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as a Gaussian distributed random variable with zero mean and known variance. We have also derived 
an approximate CRB expression when the bias PDF is very informative. Using these CRB expressions, 
we have studied the impact of the bias distribution on the MSE bound corresponding to the CRB. The 
analysis has been corroborated by numerical experiments. 
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